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The specific heat, heat flux (DTA trace) and dielectric constant of KDP ferroelectric crystal have 
been measured simultaneously for various electric fields with a conduction calorimeter. The specific 
heat presents a strong anomaly but these simultaneous measurements allow us to evaluate the latent 
heat accurately. Latent heat decreases with field and the value of critical electric field — that where 
latent heat disappears — is estimated to be (0.44 ± 0.03) kV cm -1 . Incidentally, we have measured 
simultaneously the dielectric permittivity which suggests that latent heat is developed as domains 
are growing. 
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1. INTRODUCTION 

The KDP family is one of the most extensively 
studied0,|3] hydrogen-bonded ferroelectric crystals. The 
potassium dihydrogen phosphate KH2PO4 crystal ex- 
hibits a discontinuous phase transition at To = 121 K 
from a tetragonal paraelectric phase to an orthorhombic 
ferroelectric phase. The specific heat anomaly at transi- 
tion temperature shows 0- L3- 0] a strong A- type anomaly; 
perhaps, that is why the transition was initially consid- 
ered continuous. However Reese showed that the tran- 
sition is discontinuous by measuring its latent heat, which 
was evaluated to be 46.1 J mol -1 . It is shown that dis- 
continuity disappear under the influence of pressure or 
electric field. 

The value of the critical electric field E c for which the 
discontinuity disappears is a subject of discussion. Reese 
et al.|(| also carried out a measurement of the specific 
heat with an applied field of 294 V cm -1 finding a contro- 
versial evidence of latent heat. Above E — 785 V cm" 1 
they found no evidence of latent heat and suggested 
E c = 300 V cm" 1 from the shift of the maximum of the 
specific heat as function of the electric field. Measure- 
ments of specific heat under electric field higher than 
360 V cm -1 was carried out by Sandvold and FossheimQ 
and the authors suggest that Landau theory for continu- 
ous transitions with 2 — 4 — 6 potential was appropriate 
to describe the shape of specific heat curves. 

Other works have lead to estimations of the critical 
field. For instance, Strukov et al.jg evaluated E c — 
124 V cm -1 from electrocaloric experiments. Sidncnko 
and Gladki 9] found E c — 370 V cm" 1 while Okada and 
SugiepH obtained 160 V cm" 1 for E c . Vallade[n| de- 
duced a value of 254 V cm" 1 from birefringence measure- 



ments. In contrast to these values, Kobayashi et al.[l2l] 
found, by X-ray measurements, 8000 V cm" 1 ; Eberhard 
and Horn[T3) derived a value of 6500 V cm" 1 from dielec- 
tric susceptibility. 

In this frame it would be interesting to evaluate the 
latent heat as a function of the electric field to determine 
the value of the critical electric field E c for which latent 
heat becomes null. 

The difficulty of an accurate determination of E c from 
calorimetric measurements in standard equipments - 
such as differential thermal analysis DTA and differential 
scanning calorimeter DSC — is that these systems really 
measure changes of enthalpy which has two contributions 
near the transition: one due to the latent heat and other 
due to the variation of specific heat with temperature. 
In the case of phase transitions near a tricritical point 
or a discontinuous ferroelectric phase transition under 
an electric field close to the critical field, the specific 
heat presents a strong anomaly and latent heat becomes 
very small. This fact makes difficult to separate both 
contributions and to distinguish the temperature inter- 
val where the latent heat is present. This may explain 
the lack of study of the influence of electric field, smaller 
than En, on specific heat and latent heat of KDP after 
Reese Q. 

Our group has developed a method, named square 
modulated differential thermal analysis SMDTA^M. ITU 
IT^ . based on conduction calorimetry, which is able to 
measure absolute values of specific heat and the heat 
flux exchanged by the sample when its temperature is 
changed at a rate as low as O.IKh" 1 . The comparison 
of the data allows us to separate the above two contri- 
butions to the total enthalpy and to evaluate the latent 
heat, in case there were any. 
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This technique has been successfully applied to the 
study of the almost tricritical phase transition of 
KMnF3 ^ij , whose latent heat was firstly measured with 
this technique. Furthermore, the effect of the substitu- 
tion of Mn by Ca was also investigated measuring the 
latent heat[lja], which showed that the doping makes 
the transition become continuous The method of 
SMDTA has been also applied to show that the phase 
transition in CoO[3, whose character was also contro- 
versial, is continuous. 

In this paper we have applied this method to study 
a KDP single crystal. We have measured the specific 
heat and the heat flux exchanged by the sample in the 
neighborhood of the ferroelectric phase transition. The 
measurements have been carried out at four values of E: 
V cm- 1 , 100 V cm" 1 , 400 V cm" 1 , 1000 V cm" 1 . 

Simultaneously to these measurements, the dielectric 
susceptibility of the sample has also been measured and 
we have related its behaviour around transition tempera- 
ture with the temperature interval where the latent heat 
is produced. Dielectric measurements provide informa- 
tion about the mechanism of ferroelectric phase transi- 
tion of KDP crystalfjjj [20| and the simultaneous mea- 
surement of thermal and dielectric properties would be 
worthy. 



2. EXPERIMENTAL 

The measurements were performed in a high resolution 
conduction calorimeter which has been described previ- 
ously in details pll I22I |23| . The sensor is formed by two 
identical heat fluxmeters, each one having 48 chromel- 
constantan thermocouples connected electrically in series 
but thermally in parallel. The sample is pressed between 
both fluxmeters whose signal is measured by a Keithley 
182 nanovoltmeter. Two electrodes and two heaters are 
placed between sample and fluxmeters. 

The sensor is placed inside a calorimeter block which is 
suspended within two cylindrical radiation shields. The 
whole assembly is then placed in a hermetic outer case 
at a high vacuum. The device is then surrounded by a 
coiled tube and placed in an alcohol bath. Liquid N2 cir- 
culates through the coil and regulates the temperature 
of the bath with a good thermal stabilization. As a re- 
sult, it is possible to change smoothly the temperature 
of the sample (at a rate of about 0.1 K h" 1 ) without ob- 
serving significant temperature fluctuations (always less 
than 10 -6 K) in the block temperature. 

The specific heat is measured using the method previ- 
ously described |3 ■ The same constant power W is dissi- 
pated in both heaters (dissipation branch) for twelve min- 
utes and a steady state characterized by a constant tem- 
perature difference between the sample and the calorime- 
ter block is reached. The power is then cut off until a 
new steady state is reached twelve minutes later (relax- 
ation branch). Then, the power is again switched on and 
the sequence is continuously repeated while the temper- 



ature of the assembly is changed at a low constant rate. 
That is a long-periodic serial of square thermal pulses 
is superposed to a heating or cooling ramp. The tem- 
perature increase of the sample due to the thermal pulse 
of twelve minutes (circa 50 x 10 K) is higher than the 
temperature variation of the sample produced by the rate 
of change of temperature in these twelve minutes (circa 
±20 x 10 -3 K). Hence, the sample is being cooled and 
heated alternatively during a run. 

The integration of the electromotive force given by the 
fluxmeter between every pair of steady states allows us to 
determine sample thermal capacity. Hence, the method is 
able to determine two data of heat capacity in each cycle. 
The first one is calculated from the dissipation branch 
Cd, and the second one from the relaxation branch C r . 
Heat capacity obtained in either branch show a regular 
behaviour if there is no phase transition or if it is continu- 
ous. When a discontinuous phase transition occurs, both 
data become different showing an anomalous behaviour 
in the temperature interval where the latent heat is pro- 
duced as a result of the thermal hysteresis and transition 
kinetics. That behaviour is an evidence of the discon- 
tinuous character of the transitional. Incidentally, we 
must point out that specific heat data are not reliable 
when this behaviour is observed. 

On the other hand, the DTA trace is continuously mea- 
sured in a second run without dissipation in the sample 
and using the same temperature scanning rate used to 
measure the specific heat. Due to the high number of 
thermocouples and their good thermal stability of the 
sample, the equipment works like a very sensitive DTA 
device. The electromotive force given by the fluxmeters 
is proportional to the heat flux, cf>d exchanged between 
sample and calorimeter block. 

From the specific heat data obtained in the first run- 
ning and using a method previously described|l4l[T5lll6|. 
we calculate the heat flux <p c which would have been due 
exclusively to the behaviour of the thermal capacity of 
the sample around the transition temperature. Compar- 
ing the measured tfid and the calculated <j> c we deduce that 
only in the temperature range (Tf,T p ) where both data 
do not coincide there is an effect from the latent heat. 
Its value is determined by integrating 4>d/v, where v is 
the rate of temperature change, between (Tf, T p ) and us- 
ing the straight line 4>d/v(Tf) — (f>d/v(T p ) as baseline [lfij. 
The sensitivity of the method is estimated to be better 
than 5 mJ. 

The single crystal of KH2PO4 was grown at the In- 
stitute of Physics of Poznan University (Poland). The 
sample has 0.3857 g mass, with a thickness of 2.16 mm 
along the ferroelectric axis and electrodes circular faces 
with 78.5 mm 2 in surface. The sample was placed out in 
the calorimeter. Gold electrodes were evaporated of the 
surface of the sample; those electrodes were connected to 
a capacitance bridge ESI — SP54O0 which has allowed 
us to measure the dielectric permittivity of the sample, 
simultaneously to the heat flux, with an imposed exter- 
nal bias field of V cm -1 , 100 V cm" 1 , 400 V cm" 1 and 
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1000 V cm" 1 . 



3. RESULTS 
3.1. Calorimetric measurements 

The temperature dependence of the specific heat c p of a 
sample of KDP was measured on cooling, on quasiestatic 
conditions at a scanning temperature rate of ca. v ~ 
0.1 Kh -1 using the method described in Sec. for differ- 
ent applied electric fields. In Figure ^ the specific heat 
data in a wide temperature interval for different applied 
electric field (a)£ = 0V cm" 1 , (b) E = 100 V cm" 1 , (c) 
E = 400 V cm- 1 and (d) E = 1000 V cm" 1 is shown. 

These specific heat data show a linear temperature de- 
pendence in the paraelectric phase and for V cm , 
100 V cm -1 and 400 V cm -1 a sharp A-type anomaly 
in a narrow temperature interval is observed. For 
1000 V cm -1 the maximum is more rounded and the 
transition is smeared. The shape of specific heat curve 
for high electric field tends to be almost symmetric, such 
behaviour was also suggested by Reese g|. The maxi- 
mum value of c initially increases with increasing field 
but decreases for sufficiently high values. As expected, 
the specific heat tail in the paraelectric phase increases 
with field as a consequence of the coupling of the order 
parameter to the field. 

In Figure 2, we have plotted the specific heat ex- 
cess obtained in the dissipation branch (cd) and in the 
relaxation branch (c r ) in a narrow temperature inter- 
val, 1 K, for each electric field. We observed that for 
V cm- 1 , 100 V cm- 1 and 400 V cm -1 , a and c r data 
do not coincide around the transition temperature while 
for 1000 V cm -1 both series of data almost coincide in 
the whole range of temperature. As we have stated 
above, the temperature variation of the sample due to 
the thermal pulses is slightly higher than the variation 
due to the temperature ramp. This means that in ev- 
ery period the temperature of the sample increases and 
decreases consecutively. Due to thermal hysteresis, to ki- 
netics of the phase transition etc, the process of heating 
and cooling when two phases coexist is different and con- 
sequently data obtained in the dissipation branch and re- 
laxation branch become different under these conditions. 
We have reported previously that the difference is very 
notorious even in systems near the tricritical point, where 
the latent heat is very small ( L = 0.13 J g" 1 for KM11F3, 
L = 0.010 J g" 1 for KMn . 997 Ca .oo3F 3 )[i3. 

Hence, in the case of KDP, it is clear from Figure[3]that 
for OVcm -1 , 100 V cm- 1 and 400 V cm- 1 the phase 
transition is discontinuous. On the contrary, the simi- 
lar behaviour of Cd and c r for 1000 Vcm -1 . Figure 01, 
indicate that no trace of latent heat is present, so we can 
deduce that for 1000 V cm -1 the phase transition is con- 
tinuous. Hence the critical field lies between 400 V cm -1 
and 1000 V cm" 1 . 

On the other hand, the baseline respect which the 



fluxmeter emf is integrated in order to calculate the spe- 
cific heat is the underline signal due to the temperature 
ramp imposed on the calorimeter. This signal changes 
very slowly and it does not affect to the specific heat 
measurement except when a latent heat effect happens. 
It produces a non linear variation of the baseline and 
consequently erroneous data of specific heat may appear. 
Keeping in mind that the measurements have been car- 
ried out in a cooling ramp, in Figure^ the first c r at the 
beginning of the transition which does not coincide with 
the corresponding q data is lower than the regular con- 
tribution to the specific heat. In Figure 12b we can also 
observe such a similar point but this effect is smeared. 
Finally, in Figure [2t we cannot observe any decrease in 
c r . Nevertheless, c r deviates from Cd in a small range 
of temperature as in the two previous figures. Hence, 
we can deduce that latent heat decreases with field and 
should be very small for 400 V cm -1 . 

To confirm this suggestion and to calculate the latent 
heat, we measured the DTA trace in a second run, chang- 
ing the temperature of the sample at the same constant 
rate used in the specific heat measurements to make both 
sets of data comparable. We must point out that the 
rate of change of temperature is about two orders of 
magnitude lower than the minimum value achieved in 
conventional DTA equipments. In Figure |3 we repre- 
sent the heat flux given by the fluxmeters 4>d/v (DTA 
trace) and the heat flux 4> c /v calculated from the specific 
heat data, using the method previously described, for the 
fields V cm- 1 (a), 100 V cm" 1 (b), 400 V cm" 1 (c). 

For OVcm- 1 and 100 V cm" 1 cf>d/v is higher than 
4> c /v in a very small temperature interval of about 0.05 K, 
showing the effect of latent heat. For 400 V cm _1 ^/t) is 
also higher than <fi c /v, but the difference between them 
is smaller, indicating a very small latent heat. 

The choosing of the baseline to determine the latent 
heat value for V cm -1 in graph Figure is diffi- 

cult due to the very small temperature interval where 
the latent heat is present and, consequently, the few data 
recorded inside that interval. Anyway, by considering as 
baseline the straight line between the extreme tempera- 
tures where <j)d/v and 4> c /v coincide, we obtain a value of 
43 J mob 1 . This line is also represented in Figure^). 
The value obtained for KDP deuterated at 80% using the 
same set-up and procedure was 317 J mol _1 [24|]. 

To confirm the validity of that baseline election we 
again measured the heat flux exchanged by the sample 
in identical conditions. Those heat flux data are also 
plotted in the Figure |3Ji(ii) . The figure shows that both 
heat flux data coincide with <p c (stars) in the same range 
of temperature thus supporting the baseline used for the 
determination of the latent heat. In fact, the latent heat 
for the second run is 45 J mol -1 in good agreement with 
that of the first one. On the other hand, these values 
of the latent heat are in agreement with Reese 0- It is 
noteworthy that peak area — i.e. latent heat — is quite 
reproducible despite the kinetic effects that showed each 
experiment. On the contrary, the difference observed if 
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FIG. 1: The specific heat of KDP for different appiied electric fields in a wide temperature interval. From left to right and top 



to bottom, (a) V cm _1 and (b) 100 V « 



(c) 400 V cm^ 1 and (d) 1000 V . 



(j)/v during the phase transitions suggests the existence 
of non-equilibrium kinetic process as those expected for 
a discontinuous phase transition — for instance phase 
front generation — which, as a general fact, are nonre- 
produciblc. 

For E = 100 Vcm- 1 and E = 400 V cm" 1 we de- 
creased the scanning rate of temperature when measuring 
specific heat so as to get a higher number of data points 
in the neighborhood of the phase transition. Hence, the 



determination of the baseline becomes easier as shown in 
Figure 01 The integration of the peaks gives 35 J mol -1 
and 4.2 J mol -1 respectively. 

The obtained result confirm that latent heat dimin- 
ishes with electric field (see Figure 0J). At 400 V cm -1 
latent heat is reduced by and order of magnitude with 
respect to that of zero field and it is close to the criti- 
cal field. Assuming a linear behaviour of the latent heat 
with the field we deduce that the critical field lies on 
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FIG. 2: The same as figure but in a narrower temperature interval around the transition point. Bold points stand for the 
dissipation branch (c^), open points stand for the relaxation branch (c r ). 



(0.44 ±0.3) kVcm^ 1 . 

3.2. Dielectric measurements 

In order to relate calorimetric and dielectric behaviour 
the dielectric permittivity along the ferroelectric axis has 
been measured simultaneously to the heat flux. 

In Figure we represent e(T) for the different elec- 



tric fields. The dielectric permittivity increases in all 
cases during the phase transition following the Curie law. 
In the ferroelectric phase the dielectric permittivity re- 
mains in a plateau, indicating also the large domain wall 
contribution to the dielectric permittivity, but the max- 
imum of permittivity decreases with the applied electric 
field. It has been previously reported that the maxi- 
mum value of dielectric permittivity is due to the con- 
tribution of domain wall |25|: under electric field the 
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FIG. 3: Heat flux divided by rate of change of temperature in the neighborhood of the KDP ferroelectric phase transition. 
From left to right, top to bottom, (a) V cm -1 , (b) 100 V cm -1 , (c) 400 V cm -1 . Points represent heat flux data given by the 
fluxmeters — in (a) two similar runs labelled (i) and (ii) were carried out — . Stars represent the contribution due to specific 
heat of Figure [5] Straight lines show the baseline used for determination of the latent heat which happens to be the peak area 
subtended by experimental points and the baseline. 




FIG. 4: Plot of the latent heat Ah at various electric fields. 
Critical field is found to be (0.44 ± 0.03) kV cm" 1 
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sample becomes close to monodomain state, so num- 
ber of domain walls decrease and also does the dielectric 
permittivity. Nevertheless the behaviour at transition 
temperature is different for V cm -1 and 100 V cm -1 
— discontinuous transition — than for 400 V cm -1 
close-to-critical-point transition — and 1000 V cm -1 
continuous transition — . We will relate in Figure 
the shape of dielectric permittivity for V cm^ 1 and 
400 V cm" 1 . 

The heat flux and the inverse of the dielectric permit- 
tivity versus temperature for E — V cm -1 are repre- 
sented in Figure |SJa). Three regions may be distinguish 
in this figure: (i) the paraelectric phase where the permit- 
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FIG. 5: Plot of the dielectric permittivity at various electric 
fields. Permittivity data were measured simultaneously to 
heat flux data of figure 



tivity follows the Curie law, (ii) the phase transformation 
interval where the permittivity slightly deviates from the 
previous behaviour and (iii) the ferroelectric phase where 
the permittivity shows a plateau. We must point out that 
the maximum of the permittivity matches with the end 
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FIG. 6: (a) Plot of the inverse of the dielectric permittivity 
(right axis) and ^w -1 (left axis) for E — V cm -1 , (b) Plot 
of the dielectric permittivity (right axis) and (fiv -1 (left axis) 
£ = 400 V cm" 1 . 

of phase transformation. 

For 400 V cm -1 , e(T) shows a different behaviour - 
see Figure E^b) — . At the transition temperature it 
is present a minimum. This has been observed by 
Bornarel^U in crystals of KDP for higher values of elec- 
tric field (about 1 kV cm -1 ). In that work Bornarel ex- 
plains this behaviour in terms of some domains arrange- 
ments during the phase transition; e(T) is described as 
the sum of the behaviour of £\{T) contribution that cor- 
respond to the behaviour of monodomain sample, and the 
corresponding £2(T) due to the contributions of domains. 
Moreover £i(T) increases from paraelectric phase and af- 
ter the maximum at transition temperature decreases to 
zero in ferroelectric phase; £2(T) increases from zero at 
transition temperature and remains in a plateau in ferro- 
electric phase. The sum of both contributions gives the 
appearance of e(T) in figure^. 



It may be seen also from the results presented on Fig- 
ure that the peak of the heat flux at transition tem- 
perature appear at temperature where the production of 
domains become dominant. 



4. CONCLUSIONS 

The thermal and dielectric behaviour of KDP crys- 
tal near the temperature of its ferro-paraelectric phase 
transition has been simultaneously studied under the in- 
fluence of electric field. For V cm -1 the latent heat 
has been measured twice. Although both measurements 
shows a different kinetic, their results are similar and in 
good agreement with those obtained by ReeseQ. Despite 
the high increase of the specific heat around the tran- 
sition temperature, the very small values of the latent 
heat and the narrow temperature range — ca. 0.1 K — 
where transition is developed, we have been able to dis- 
criminate the contribution of the latent heat to the total 
change of enthalpy. Finally, simultaneous measurement 
of heat flux and dielectric susceptibility suggest that the 
effect of latent heat appears at temperatures range where 
domains are growing. Maximum dielectric permittivity 
was noticed at temperature where phase transition ends. 

It was established that the latent heat decreases with 
the field and the critical electric field is estimated to be 
(0.44 ±0.03) kV cm -1 . 
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